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Abstract 

The low energy theorems for QCD are generalized to finite tem- 
perature and chemical potential, including non-zero quark masses. 
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Low-energy theorems for quantum chromodynamics (QCD) at zero tem- 
perature and density were derived long ago by Novikov et al. [Ij]. They are 
useful in a number of contexts, for instance in constraining effective theo- 
ries or in assessing lattice gauge calculations. Recently the theorems were 
generalized to finite temperature (T > 0) for the pure glue sector of QCD 
[0]. The purpose of the present note is to provide a further generalization 
by including the quark sector and allowing for non-zero density or chemical 
potential, /i. 

For clarity we shall consider one quark flavor since the case of several 
flavors is trivially obtained by introducing flavor-dependent masses, mo, and 
chemical potentials, /i, and appropriately summing over the flavor index. In 
the imaginary time approach the partition function takes the familiar form 
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[dA][dq][dq] 
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Here we have suppressed the gauge fixing and Fadeev-Popov ghost terms, as 
well as the quark color labels since they are inessential here. The generators 
of color SU(3) are denoted by A a , and the gluon fields and field strength 
tensors have been scaled by the bare coupling constant g : A£ = QqA 1 ^ and 

F % = 9oF%v 

We first consider the case where the quark mass is set to zero. The grand 
potential of the system is defined in the usual way, Q = —T In Z, and we 
have 

& 9o /rpuu/n n \ pa lr\ n\\ — 9o 



(*T (0, 0)^(0, 0)) = (F 2 ) , (2) 



d(-l/gi) V 4 \ a K ' ' » uK ' ; / 4 

where V is the volume of the system. The angle brackets denote a thermal 
average. 

This derivative can be calculated in another way by using an explicit form 
for Q/V determined on dimensional grounds. Regulation of the ultraviolet di- 
vergences of the theory introduces a mass scale, M, at which the subtractions 
are performed. The results can be written in terms of the non-perturbative 
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dimensionful parameter 



A = Me W { I > ( 3 ) 



x /3 s (oc s 

(a s (M) 

where a s = <?o/47r and f3 s is the Gell-Mann-Low function: Mda s /dM = j3 s . 
There are two additional dimensionful parameters, namely T and /i. Since 
the grand potential is an observable quantity with zero anomalous dimension 
[[J, we can write in general 

where / is some function. We note that at zero temperature and chemical 
potential a form proportional to A 4 can be formally justified within a well- 
defined regularization scheme 0. In Eq. (|J) g is involved only through A, 
hence we obtain 

d a 4™? / d d\n 

4 - T 7^-/^7 T7- ( 5 ) 



d(-l/g 2 )V (3 s (a s )\ dT ^dfiJV 
This leads to the chain of equations 

We have used the standard QCD result @] for the trace of the improved 
energy-momentum tensor density 6^ |J. We have also used the standard 
thermodynamic relation for the energy density 

and the pressure P = —Q/V. 

We turn now to the case where the quark mass is non-zero. Since Q is a 
physical observable it has to be expressed in terms of renormalization group 
invariant quantities. However, the quark mass has anomalous dimension and 
depends on the scale M. The renormalization group equation for mo(M), 
the running mass, is (M / 'm^jdrriQ / 'dM = — 7 m and we use the MS scheme for 
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which (3 S and 7 m are independent of the quark mass ||. Upon integration 
the renormalization group invariant mass is given by 



a s {M) 



rh = itlq{M) exp 



7m (« S ) 



da* 



where the indefinite integral is evaluated at a s (M). Then Eq. 

A 4 h 
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becomes 
(9) 



where h is some function. Proceeding as before we obtain 
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Here we have used the trace of the energy-momentum tensor for QCD with 
quarks 0, Clearly the physical quantities in Eq. flltf) must obey the same 
relation as in Eq. (§). 

We can iterate this procedure by taking n further derivatives of Eq. fllOD . 
In doing so it is convenient to note that since the grand potential density 
and its derivatives are independent of the scale M at which the ultraviolet 
divergence are regulated, we can choose any scale to prove the result. It is 
convenient to pick a sufficiently large scale that we can take the lowest order 
expressions, (3 s {a s ) -> -ba 2 s /2ir, where b = (UN - 2N f )/3 for the SU(N) 
gauge group with Nf flavors, and (1 + j m ) — > 1. Then it is straightforward 
to obtain the following relation 



n+l 
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dr n d 6 x n ■■■ dr 1 d i x l (9£(T n , Xn) • • • 0£(ri, xi)^(0, 0; 



Here the subscript c means that only the connected diagrams are to be in- 
cluded and the limits of the imaginary time integrations are suppressed. It is 
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possible to perform a similar analysis for a renormalization group invariant 
operator O 

( T ^ + 4- d )" <0> 

= J dr n d 3 x n ■■■ J dnd 3 Xl (0£(r n , x n ) • • • 0£(n, Xi)0(O, 0)) c , (12) 

where d is the canonical mass dimension of O. If the operator O also has 
anomalous dimension then the corresponding 7-function will have to be in- 
cluded. Equations fllTD and (|T^) differ from those found previously [Q] by 
the addition of the operator fid/dfi. One could also similarly generalize the 
finite-momentum low-energy theorems [Q, [J. 

As an example, where we can evaluate the left hand side of Eq. ([TT|), 
consider the case where /1 and/or T is much greater than A so that per- 
turbation theory is valid. Considering the first term in the /^-function, the 
behavior of the strong coupling constant as a function of chemical potential 
and temperature can be written 

a - (f) = bim • (13) 

where £ = Ty(fi/T) and we need not specify the function y. The perturbative 
expression for the pressure in SU(N) gauge theory up to two-loop order is 



P = —(N 2 - 1)T 4 + N —— + + 



7T 2 _ ^ ,./7vr 2 T 4 , /i 2 T 2 , /x 4 



45 v ~ _/ ' ~ V 180 ' 6 ' 12tt 2 i 

<N2 ~ 1] a.(t) ( (4iV + 5)T 4 + !L 2 T 2 + ^ ' ) . ( I I. ) 



Then from Eqs. flT3|) and ( |14"D we find 



<9T ^«9/i / \ "/ 288 
x<+ 2 (0 (^V(n + 1)! ((4iV + 5)T 4 + ^ 2 T 2 + A//) . (15) 

This provides a requirement on model or lattice calculations of the right hand 



side of Eq. (11). Note, however, that if we consider the next order in the 
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/3 s -function then, for N = 3 and one flavor in the MS-scheme, we obtain a 
correction factor of order (1 + lAa s n) to Eq. flTop . Thus, at this level of 
approximation for P and a s , an arbitarily large number of derivatives cannot 
be taken without at some point losing all accuracy. 

In conclusion we have shown that at finite temperature and chemical 
potential the low-energy theorems of Novikov et al. hold provided that 
the operators Td/dT and fid/dfi are appropriately included as in Eqs. (|TT|) 
and (0). 
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